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00 : ABSTRACT 

We report on the study of a sample of 25 field early-type galaxies, in the redshift 
' range z « 0.1 — 0.5, selected on the basis of colours and morphology from the Medium 

Deep Survey (MDS) of the Hubble Space Telescope (HST). Surface photometry in two 
colors (F606W and F814W) and redshifts have been derived for all the galaxies in the 
sample, while velocity dispersions has been measured for 19 of the sample galaxies, as 
described in a companion paper. 

Our study of the evolution of the Fundamental Plane (FP) with redshift yields 
, the following results. Field early-type galaxies define a tight FP out io z ~ 0.4, with 

' scatter unchanged with respect to local samples, within the observational errors. The 

' intermediate redshift FP is offset with respect to the local FP of the Coma Cluster, in 

the sense that, for given effective radius and velocity dispersion, galaxies are brighter 
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Q ■ than expected from the local relation. The offset of the FP is found to increase with 

^ \ redshift. The range of parameters covered by our sample is not sufficiently extended 

C/3 i to measure the slopes of the FP at intermediate redshift. Similar results are found for 

^ ' the SBq-Rq relation, out to z « 0.5. 

The evolution of the FP (and of the SBq-Rq relation) is studied quantitatively 
with a Bayesian-Montecarlo technique. By applying this technique, we find that the 
offset of the intercept of the FP (A7) with respect to the local FP increases as A7 = tz 
with the following 68 per cent limits: 0.33 < r < 0.44 (for = l,i^A = 0) or 
0.44 < T < 0.56 (for Vl — 0.3, J^a = 0.7). In addition, we interpret the results in 
terms of the evolution of the stellar populations, under the assumption of passive 
evolution. In a single-burst scenario, the observed properties are consistent with those 
of a stellar population formed at 2; ^ 2 (for = 1, JIa = 0, Hq = 50 km s~^Mpc~^) 
or 0.8 < z < 1.6 (for = 0.3, I^a = 0.7, = 65 km s^^Mpc"!). If a small fraction of 
the stellar mass is formed in a secondary burst, the primary burst may have occurred 
at higher z. For example (for Q, ~ 0.3, J^a — 0.7, Hq — 65 km s^^Mpc^^) the primary 
burst may have occurred at z ^ 3 if a secondary burst with a tenth of stellar mass 
occurred at 0.6 < z < 0.8. 

Finally, the intercept and scatter of the FP found for field early-type galaxies and 
for cluster data (taken from the literature) at z « 0.3 — 0.4 are mutually consistent, 
within the observational errors. If higher redshift (up to z = 0.83) cluster data are 
considered, the ages of the stellar populations of field early-type galaxies inferred from 
a single-burst scenario are found to be marginally smaller than the ages derived for 
the cluster galaxies. 

Key words: 
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1 INTRODUCTION 

Despite the progress made in understanding the physics of 
early- type galaxies (E/SO), a single widely accepted model, 
capable of explaining all their observed properties, is not 
yet available. A variety of scenarios are still considered to 
be viable (see e.g. Larson 1975, Toomre 1977, van Albada 
1982). Schematically, we can identify the following two main 
classes of scenarios. In the traditional view, the monolithic 
collapse scenario, E/SO form in a rapid collapse shortly after 
an intense burst of star formation. Since the stellar popula- 
tions of E/SO in the local Universe are generally old, and 
no such intense bursts of star formation are observed at 
redshift z < 1, it is generally assumed that most of these 
monolithic collapse events occurred at high redshifts. In the 
alternate view, the merging scenario, E/SO form by interac- 
tions of disc galaxies. The merging scenario has been cast 
into a cold dark matter cosmological framework, the hier- 
archical clustering scenario (White & Rees 1978; Kauffman, 
White & Guiderdoni 1993). In the hierarchical clustering 
scenario the formation of massive E/SO by merging of discs 
is a relatively recent phenomenon (Kauffman 1996). 

In the absence of a satisfactory theoretical scenario, new 
observations play a key role, by providing additional input 
and the opportunity to test predictions. Since the theoreti- 
cal models have been conceived and 'calibrated' in order to 
reproduce the properties of local galaxies, the observation 
of intermediate (0.1 < z < 1) and high redshift {z > 1) 
galaxies is particularly important. 

Several observational results have been obtained on this 
subject in recent years. Most of them have been obtained for 
the cluster environment, where early- type galaxies are more 
abundant and studied more efficiently (Ellis et al. 1997; 
Dressier et al. 1997; Stanford, Eisenhardt & Dickinson, 1998; 
van Dokkum et al. 1998a,b; Brown et al. 2000; Kelson et al. 
2000b; van Dokkum et al. 2000) . These studies indicate that 
the stellar populations of early-type galaxies in the core of 
rich clusters have possibly formed at high redshift, earlier 
than z ~ 2 — 3. 

When comparing intermediate redshift E/SO to local 
ones, one should take into account the evolution with red- 
shift of clusters and of the population of galaxies in clusters. 
Although the evolution of clusters of galaxies is still some- 
what controversial (Postman et al. 1996), it is is likely that 
the population of galaxies in clusters evolves by encroach- 
ing on isolated galaxies and small groups. Thus, part of the 
progenitors of present-day cluster E/SO were likely located 
outside clusters in the past. Therefore, to gain a reliable 
global picture, it is necessary to study the evolution of E/SO 
in all environments, from the core of rich clusters to the 
field. In addition, E/SO have been observed to be subject to 
interactions with their environment ranging from weak dis- 
tant encounters to major merging events. Thus, since the 
frequency of interactions depends upon the environment, 
evolution should depend upon the environment as well, if 
interactions play a significant role. For this reason, the role 
of environment should be especially important in the merg- 
ing scenario. In the hierarchical clustering scenario, numer- 
ical simulations predict (Kauffmann 1996) that field E/SO 
should have younger stars than cluster ones and have as- 
sembled later. Therefore, the comparison of the evolution of 



the properties of cluster and field E/SO provides a test of a 
prediction of hierarchical clustering models. 

So far, only a few studies of field early-type galaxies, 
based on HST high resolution images, have been published 
(Im et al. 1996; Treu & Stiavelh 1999; Schade et al. 1999; 
Menanteau et al. 1999; Brinchmann & Ellis 2000), often with 
very limited statistics (e.g. the analysis of early-type galaxies 
in the Hubble Deep Fields, Franceschini et al. 1998; Benitez 
et al. 1999; Kodama, Bower & Bell, 1999; see also Treu et al. 
1999). One important result is that field early- type galaxies 
at intermediate redshifts span a wider range of colours than 
cluster ones (more extended to the blue; Schade et al. 1999). 
This is likely to result from a more complex evolution of the 
stellar populations, with episodes of star formation occur- 
ring at lower redshifts than in clusters (see e.g. Trager et al. 
2000a,b). A detailed study of the star formation history of 
field E/SO is needed to clarify this issue. 

A correct picture of the evolution of stellar populations 
is also needed in order to understand how the number den- 
sity of E/SO evolves with redshift. In fact, several studies 
(Treu & Stiavelli 1999; Menanteau et al. 1999) find that the 
number of red morphologically selected E/SO at high z is 
smaller than what is predicted by models based on passive 
luminosity evolution and constant comoving density mod- 
els. However, Schade et al. (1999) find that the number 
density of early-type galaxies, selected purely on morphol- 
ogy without color criteria, is consistent with being constant 
from 2 ~ 1 to 2 ~ 0.2 (see also Im et al. 1996). Similarly, 
Menanteau et al. (1999) find that, if all E/SO are considered 
(including the blue ones), the data are consistent with no 
evolution of the number density. 

The evolutionary history of E/SO derived from the data 
depends critically on the detailed star formation history in- 
volved (for other sources of uncertainties see, e.g., Chiba & 
Yoshii 1999; Schade et al. 1999; Barger et al. 1999; Daddi 
et al. 2000). In fact, even otherwise negligible amounts of 
young stars are sufficient to make high— 2; E/SO appear blue 
and thus to produce the observed "deficit" of red high-z 
E/SO (Treu & Stiavelh 1999; see also Jimenez et al. 1999 
and Abraham et al. 1999). 

More detailed spectroscopic information is needed to 
study the evolution with redshift of the structural proper- 
ties of E/SO and of their stellar populations. For this reason 
we have initiated (Treu et al. 1999; hereafter T99) an ob- 
servational project aimed at obtaining high signal-to-noise 
medium resolution spectra for a sample of field early-type 
galaxies at intermediate redshift selected on the basis of 
HST-WFPC2 images. The photometric and kinematic mea- 
surements for the entire sample (0.1 ^ z 0.6) are described 
in a companion paper (Treu et al. 2001; hereafter PII). In 
the present paper, the third of the series, we study the evo- 
lution with redshift of the Fundamental Plane (Dressier et 
al. 1987; Djorgovski and Davis 1987), a very tight empiri- 
cal correlation between the main observables of E/SO (Sec- 
tion ^). Its existence and tightness have profound implica- 
tions in terms of galactic formation and evolution, and it is 
therefore particularly interesting to investigate how far in 
the past it exists and what is the time evolution of its co- 
efficients and scatter. In addition, the Fundamental Plane 
is used as a diagnostic of stellar populations, to study the 
star formation history of field E/SO. A similar study is per- 
formed on the evolution of the SBe-Re relation (Kormendy 
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1977; Section |4.2| ) , a photometric projection of the Funda- 
mental Plane. Future papers will report on the measurement 
of ages and metallicity, using the absorption line indices in 
the Lick/IDS system (Trager et al. 1998). 

The paper is organized as follows. The Fundamental 
Plane as a diagnostic tool of the stellar populations is de- 
scribed in Section ^. In Section ^ we give a brief summary 
of the data described in PIl. In Section ^ we illustrate the 
observational results, i.e. we discuss the existence and loca- 
tion of the intermediate redshift field Fundamental Plane. 
In Section ^ we analyse selection effects that must be taken 
into account while interpreting the results in terms of stellar 
population evolution (Section P). In Section ^ we compare 
the Fundamental Plane of field early-type galaxies to the 
Fundamental Plane of clusters at similar redshift, with data 
and results taken from the literature. Conclusions are drawn 
in Section ^ 

The results presented here are given for two sets 
of values of the cosmological parameters (the matter 
density and the cosmological constant SIa, expressed 
in dimensionless form, and the Hubble constant Ho = 
50/i5okm s~^Mpc~^): a long lived A Universe (f2 — 0.3, 

= 0.7, /i5o=1.30; age of the Universe ~ 14.6 Gyr) and 
the 'classical choice' of parameters {Q = 1, JIa = 0, /150 = 1; 
age of the Universe ~ 13.1 Gyr). We will refer to them as 
the A cosmology and the classical cosmology respectively. 



2 THE FUNDAMENTAL PLANE OF 
EARLY- TYPE GALAXIES 

Early-type galaxies are a rather homogeneous family of 
galax;ies. It has been early-on recognized that several ob- 
servables, describing their dynamical state, size, and chemi- 
cal composition, correlate with luminosity. Typical examples 
are the effective radius (i?e; Fish 1964), the central velocity 
dispersion (ct; Faber & Jackson 1976), the effective surface 
brightness {SBe; Binggeli, Sandage & Tarenghi 1984), the 
integrated colours (Bower, Lucey & Ellis 1992 and references 
therein), and metal indices such as Mg2 (Bender, Burstein 
& Faber 1993 and references therein). 

The relatively large scatter of these correlations sug- 
gests that a 'second parameter' could be involved. The 'sec- 
ond parameter' was found when large spectrophotometric 
surveys during the mid-eighties revealed the existence of 
a tight relation between effective radius, effective surface 
brightness, and central velocity dispersion (Djorgovski & 
Davis 1987; Dressier et al. 1987). The relation, 

log iie = Q log a + (3 SBe + 7, (1) 

where Re is in kpc (we will reserve the notation re when 
such radius is measured in arcsec), a in km s^^, and SBe 
in mag arcsec"'^, was called the Fundamental Plane (FP). 
Because the value of the Hubble constant enters the relation 
through the calculation of the effective radius in kpc, the 
value of 7 depends on Ho- In the following we will call FP- 
space the three-dimensional space defined by logi?e, logo-, 
and SBe. 

The quantities SBe and Re are found to depend on the 
wavelength of observation, as a result of colours and colour 
gradients in the stellar populations of early-type galaxies. 
When needed, we will explicitly state the dependence on 



wavelength by adding a subscript (e.g., -Res is the effective 
radius in the B band). Not much is known about systematic 
variations of a with the spectral range used for the mea- 
surements. In PII we showed that no significant variation is 
noticed if lines such as Ga H and K or NaD are included 
or excluded from the spectral region of the measurement. 
In the following, in order to compare our data to data from 
the literature, we will assume that the value of the central 
velocity dispersion does not depend on the spectral region 
used for the measurement. 

The intercept 7 depends on wavelength. For example, 
in the Goma Gluster, Bender et al. (1998) report in the B 
band 7s = —8.895 assuming /150 = 1; in the V band, using 
the sample and slopes of Lucey et al. (1991), 7v = —8.71 
is found. The determination of 7 is hampered by the errors 
in the determination of distances and extinction. No dif- 
ferences, within the observational errors, have been found 
among the intercepts (and slopes) of the FP of nearby clus- 
ters and the field (Pahre, de Carvalho & Djorgovski 1998). 

The slope /3 is very well determined, independently of 
the different selection criteria and fitting techniques adopted 
by various groups; in addition, it does not change signif- 
icantly with wavelength (Pahre, Djorgovski & de Garvalho 
1998; Scodeggio et al. 1998). A typical value is /3 = 0.32. The 
slope a suffers from larger uncertainties and it is observed to 
vary with wavelength (Pahre et al. 1998b ; Scodeggio et al. 
1998). Accurate measurements over large samples of galaxies 
have shown that the observed scatter of the FP is larger than 
the observational errors (e.g. j0rgensen, Franx & Kjargaard 
1996). In other words, a significant fraction (~ 0.07 — 0.08 in 
log _Re) of the observed scatter is due to scatter in the phys- 
ical properties of E/SO. Within the observational errors, the 
intrinsic scatter is observed to be constant with wavelength 
and along the FP. 

2.1 The physical origin of the FP 

The very existence of the FP is a remarkable fact. Any the- 
ory of galaxy formation and evolution must be able to ac- 
count for its tightness and its universality. In the following 
we will address the issues of how the observables entering 
the FP may be related to intrinsic physical quantities such 
as mass and luminosity of early- type galaxies and we will 
list some of the physical processes that have been proposed 
to be at the basis of the FP. In addition we will describe 
how the existence of the FP and its tightness can be used to 
derive information on the formation and evolution of early- 
type galaxies. 

Let us define an effective mass M, 

M = a^Re/G, (2) 

where G is the gravitational constant. If homology holds, i.e. 
early-type galaxies are structurally similar, the total mass 
Ai (including dark matter if present) is proportional to M 

M = CiM (3) 

and we can interpret Eq. ^ in terms of the Virial Theorem. 
If ci depends on mass, then we will say that weak homology 
holds. 

Let us also define the effective luminosity L as 
- 2.5 log L = SBe - 5 log Rc-2.5 log 2tt. (4) 
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In general, from the FP and Equations Q and |4| we have: 



w — 



(5) 



The observed values for 10/3 — a — 2 cluster around to zero, 
although there is significant scatter that could be related to 
differences in the selection of the samples, in the determi- 
nation of distances, in the fitting techniques, and extinction 
corrections. If 

10/3 - Q - 2 = 0, (6) 

then Equation ^ can be reduced to a power law relation 
between effective mass and luminosity (e.g. Faber et al. 1987; 
van Albada, Berlin & Stiavelli 1995): 

L (X M\ (7) 
where 



?7 = 0.2 



(8) 



Based on these definitions, the M/L (effective mass-to-light 
ratio) of a galaxy is readily obtained in terms of the FP 
observables: 



M/L oc 10° *^^e^^^e 



(9) 



By comparing M/L to the value predicted by the FP for the 
same effective radius and velocity dispersion, 

( M/L)fp oclO"Wcr 50 R^, (10) 



the scatter of the FP can be connected to the scatter in M/L 
of the early-type galaxy population: 

\ 2 



E. [log(M/L') - log(Af/L)', 



7 



■7 



(11) 



(2.5/3)2 

where the superscript i labels the galaxies of the sample 
under consideration, and 7' is defined as 



7 



log Re — a log a' — PSBlg 



(12) 



Thus, the remarkable tightness of the FP implies a very low 
scatter in the M/L of early-type galaxies. 

When a theoretical effort to understand this empirical 
scaling law is made, it is necessary to model the relation 
between effective mass and total mass {M) and the relation 
between stellar (M*) and total mass. For example, if we 
assume 



Moi M oc M,, 

then Equations | 

M* _ T ^ 

-— cx 10 =f* Rp 



I and ^ yield 



10^- 



(13) 



(14) 



By neglecting the mild dependence on effective radius. Equa- 
tion h4j implies that the stellar mass-to-light ratio depends 
on the total mass of the galaxy. This phenomenon is some- 
times referred to as the tilt of the FP (Renzini & Ciotti 
1993; see also Bender, Burstein & Faber 1992; Pahre et al. 
1998b). 

Several interpretations in terms of properties of the stel- 
lar populations have been suggested to explain the tilt of the 
FP (see, e.g., Renzini & Ciotti 1993; Ciotti, Lanzoni & Ren- 
zini 1996; Graham & CoUess 1997; Pahre et al. 1998b). How- 
ever no conclusive interpretation of the FP has been found 



so far. Some of the hypotheses have turned out to be insuf- 
ficient to explain the FP (e.g. the tilt cannot be attributed 
entirely to a relation between metallicity and mass, Pahre et 
al. 1998b). On the other hand, the viable hypotheses require 
an 'unnatural' fine tuning in order to explain the tilt while 
preserving the tightness. 

The remarkable tightness of the FP can be used to 
set constraints on the formation and evolution of early-type 
galaxies. For example, assuming weak homology and a pro- 
portionality between total mass and stellar mass, we know 
from the FP that at any given mass the scatter in log(A/»/L) 
is very small (~0.10). If we attribute such scatter entirely to 
scatter in age, we can set an upper limit to it (At/t). Simi- 
larly, we can set an upper limit to the scatter in metallicity 
(AZ/Z). Unfortunately, correlations between the various ef- 
fects can conspire to keep the scatter small, so that, for 
example, a small scatter in \og{M^, / L) does not necessarely 
imply a small scatter in both age and metallicity. 

The absence of a difference between the FP of field and 
cluster E/SO is also a remarkable fact. For example, within 
the scheme described in this section, it can be interpreted in 
terms of the absence of significant differences between the 
average M/L, and thus between the stellar populations, of 
field and cluster E/SO. This fact is consistent with the very 
small difference found between the Mg-cr relation (in which 
the observables are distance independent) of a large sample 
of local ellipticals in the cluster and in the field environments 
by Bernardi et al. (1998). 

A very effective way to further investigate the origin 
of the FP, and thus to learn about the evolution of E/SO, 
is to study the FP as a function of redshift. In fact, with 
this information it is possible to distinguish between models 
that predict an evolution of the slopes with redshift, e.g. the 
model where the tilt is due to Initial Mass Function (IMF) 
effects (see Renzini & Ciotti 1993), and models that do not, 
e.g. the model where the tilt is due to dark matter distribu- 
tion. Similarly for the scatter: if the scatter is due to scatter 
in the ages of the stellar populations, then it should increase 
with look-back time; if it is due to scatter in metallicity, it 
should not. 



2.2 The FP at intermediate redshift as a 
diagnostic of the evolution of stellar 
populations 

Recently it has been observed that early-type galaxies in 
clusters at intermediate redshift also populate a tight FP 
(van Dokkum & Franx 1996; Kelson et al. 1997; van Dokkum 
et al. 1998b; Bender et al. 1998; Pahre 1998; j0rgensen et 
al. 1999; Kelson et al. 2000b). 

The intercept of the FP (7) is observed to vary as a func- 
tion of redshift. So far, no dramatic change in the slopes with 
redshift has been reported, even though an accurate mea- 
surement is difficult, because of the small samples available 
and of the intrinsic problems in this measurement, such as 
completeness corrections and fitting technique (see, e.g., van 
Dokkum & Franx 1996; T99; j0rgensen et al. 1999, Kelson 
et al. 2000b). For example. Kelson et al. (2000b), using the 
largest sample at intermediate redshift (2 — 0.33) published 
at the moment of this writing (30 E/SO galaxies in a single 
cluster), find a = 1.31 ± 0.13 and (3 = 0.344 ± 0.10 in the 
V band; these values are indistinguishable from those found 
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locally. The absence of a dramatic change in the slopes ar- 
gues against the interpretation of the FP solely in terms of 
an age-mass relation. 

A simple explanation of these findings can be given by 
recalling the correspondence of the FP observables to A/ 
and M/L. Let us consider a sample of galaxies with observ- 
ables labelled by superscript i. Let us assume that effective 
radius and central velocity dispersion do not change on the 
observed time-scale. For each galaxy, we can compute the 
variation of the logarithm of the effective mass-to-light ratio 
with redshift, 

A log(M/L)' = log(Af /L):=,, ~ log(Af/L)U.i , (15) 

in terms of the variation of combinations of the coefficients 
of the FP with z from Eq. |l^: 

Alog(A//Lr= A (124^) log 

+A(^)logRe-A(^), (16) 

where the symbol A indicates the difference of the quantity 
at two redshifts as in Equation 

For the analysis presented here we will assume that a 
and P are constant. This assumption is broadly consistent 
with the observations and makes the interpretation of the 
results straightforward. In fact, if a and /3 are constant, we 
have that: 



Alog(M/r)=-|g, 



(17) 



i.e. the evolution of log(A'//L*) depends only on the evolu- 
tion of 7*. In practice, by measuring 7' for a sample of Ng 
galaxies at intermediate redshift, and by comparing it to the 
value of the intercept found in the local Universe (70), we 
can measure the average evolution of \og{M/L) with cosmic 
time: 



< A log(Af/L) >= 



i=JV„ 7 



■ 70 



2.5PNg ' 



(18) 



i.e., if the intercept 7 is computed as the average of 7*, 
A7 



< Alog(A//L) - 



2.5/3' 



(19) 



Similarly, the intrinsic scatter of AI/L is related to observ- 
able quantities by Equation |ll], and therefore its evolution 
with redshift can be measured. Note that, by assuming con- 
stant Re and a, we are neglecting any structural and dy- 
namical evolution of the galaxies; in addition, by assuming 
constant a and /3, we obtain that the evolution of M/L is 
independent of M. 

In the spirit of this observational study, we will focus 
our attention on the measurement of the evolution of the 
effective mass-to-light ratio M/L and we will not attempt to 
address the issues associated with the relation between the 
effective mass M and the total mass A4. In particular, we 
will limit our study to the evolution of stellar populations, 
neglecting any dynamical or structural evolution. In other 
words, the only properties of the galaxies that will be allowed 
to vary in our study will be those of the integrated stellar 
population. The comparison with the expectations of passive 
evolution models will be done by assuming that: 



Two effects must be kept in mind when applying this 
diagnostic: 

(i) Morphological evolution. If galaxy types evolve 
into one another, then local early-type galaxies are not sim- 
ply the evolution of intermediate redshift ones. For example, 
if spirals evolve into early-type galaxies as suggested by hi- 
erarchical clustering scenarios, then the FP at intermediate 
redshift can look tight and slowly evolving because E/SO 
are being recognized only when their stellar populations are 
already old. 

(ii) Cluster population evolution. The population of 
galaxies within clusters is likely to evolve with redshift. If 
this happens, then the comparison of cluster and field sam- 
ples can be biased. 

A complete analysis of these effects requires a compre- 
hensive theoretical model, able to describe the joint evolu- 
tion of galaxies and of their environment. Such modeling is 
beyond the aims of this observational study. Therefore, we 
will mostly limit our presentation to the observational re- 
sults (to be used as a constraint when testing models). In 
order to facilitate comparison with theoretical models, the 
sample selection process has been defined as clearly as pos- 
sible in PII. Some considerations on the population biases 
in the context of simple passive evolution of the stellar pop- 
ulations are presented in Section 3.5. 



The comparison of the stellar populations in cluster and 
field environment is particularly interesting. In fact, as men- 
tioned in the Introduction, one of the expectations of hier- 
archical clustering models is that early-type galaxies and 
their stellar populations form much later in the field than 
in the core of rich clusters. Unfortunately, in the local Uni- 
verse, in addition to the observational uncertainties listed 
in Section ^, the detection of a possible difference between 
cluster and field galaxies with the FP is hampered by the 
generally old age of their stellar populations. It is easier 
to detect any possible difference at intermediate redshift, 
where the FP is in principle more sensitive. In fact, since 
peculiar velocities (the motion with respect to the local co- 
moving system) are limited (a few hundreds of km s~^, see 
e. g. Groth, Juszkiewicz & Ostriker 1989), the uncertainty 
on distance determination becomes negligibly small with in- 
creasing redshift. In addition, if the stellar populations of 
field and cluster E/SO had different average ages, any differ- 
ence in the FP should increase with redshift thus becoming 
detectable. 



2.3 The SBe-Re relation 

The correlation associated with the projection of the FP on 
the photometric plane, also known as the SBe-Re relation 
(Kormendy 1977), 



SBe = Hi log_Re + ci2 



(21) 



M oc M,. 



(20) 



is very useful for our purposes because it does not involve 
internal kinematics. Therefore, it can be followed to higher 
redshifts and fainter absolute magnitudes. The drawback is 
that the SBe-Re relation has a larger intrinsic scatter than 
the FP and its coefficients are much more sensitive to sam- 
ple selection biases (see, e.g., Capaccioli, Caon & D'Onofrio 
1992). In the following we will use the evolution of the SBe- 



© 0000 RAS, MNRAS 000, 000-000 



6 T. Treu et al. 



Re relation to extend the findings achieved with the FP to 
higher redshift and to fainter absolute luminosities. 



3 THE DATA 

The galaxies have been selected on the basis of HST-WFPC2 
images taken from the HST- Medium Deep Survey (Griffiths 
et al. 1994), mainly based on colour, magnitude, and mor- 
phology. The sample selection criteria are extensively dis- 
cussed in PII and are taken into account in the analysis 
presented here. Briefly, we can say that the sample analyzed 
in this paper is representative of the population of early-type 
galaxies at z « 0.1 — 0.5 with 0.95 < V606-l8i4 < 1.7 and 
I8i4 < 19.3 (as in PII Veoe and I8i4 indicate Vega magni- 
tudes through HST fihers F606W and F814W respectively). 

The sample is a field sample in the sense that the galax- 
ies are selected on random WFPC2 pointings, excluding the 
regions of known clusters. On the basis of the images avail- 
able, it is not possible to identify and to exclude members 
of groups or galaxies in the outskirts of rich clusters. In this 
respect, our operational definition of field does not coincide 
with the definition used in theoretical investigations, where 
the full three dimensional and dynamical information avail- 
able is used to identify rich clusters, groups, and field (e.g., 
Kauffmann 1996). 

The data, the data reduction, and the caveats of the 
photometric and spectroscopic measurements are described 
in PII. For quick reference, we summarise in Table ^ the 
relevant observed quantities for the sample of E/SO analyzed 
in this paper. 



4 OBSERVATIONAL RESULTS 

In this section we will discuss the location of our sample 
of field early-type galaxies at intermediate redshift in the 
FP-space in comparison to the location of local samples. 
Only galaxies that obey the selection criteria listed in PII 
are considered. 

In principle, one would like to compare the intermedi- 
ate redshift data to field local data. However, as discussed 
in Section |^, no significant difference has been found so far 
in the local Universe between the FP in the field and that 
in the cluster environment, within the accuracy allowed by 
the uncertainties in the determination of distances, redden- 
i ng corrections, and the small size of the available samples 
(Pahre et al. 1998a). For these reasons, we will compare the 
intermediate redshift data to the FP of the Coma cluster, a 
well defined relation, with selection biases under control. 

The discussio n wi ll be extended to a larger sample of 
galaxies in Section 



4.2 



by means of the projection of the FP 
on the photometric plane, the SBe-Re relation, including 
galaxies without measured velocity dispersion. 



4.1 The Fundamental Plane of field early-type 
galaxies at intermediate redshift 

In Figures |l| to ^ panels (a), (b), and (c), we plot the location 
in the FP-space of the intermediate redshift galaxies binned 
in redshift. The data are shown in the rest-frame B and 
V bands and for two choices of values for the cosmological 
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Figure 1. Fundamental Plane in the rest-frame B band at inter- 
mediate redshift. In panels (a), (b), and (c) we plot the interme- 
diate redshift galaxies (pentagons) binned in redshift. The solid 
line represents the FP of Coma, as measured by Bender et al. 
(1998). In panel (d) the average offset of the intercept 7 (at fixed 
slopes a and /3, taken from Bender et al. 1998) is plotted as a 
function of redshift. The error bars are the standard deviation of 
the mean. The classical cosmology is assumed. 

parameters (see captions). Since errors on the photometric 
parameters are correlated, we do not show them in this pro- 
jection. A better visualization of the errors is given in the 
projection that separates kinematics and photometry, used 
in Figures ^ and ^. In panels (d) of Figures ^ to ^ we show 
the evolution of the intercept 7 of the FP with redshift. 
Qualitatively, two main facts emerge: 

(i) Existence of the FP. At any given redshift between 
z — and z ~ 0.4 the FP of field early-type galaxies is well 
defined. The scatter is small. No trend of increasing scatter 
with redshift is noticeable within the accuracy allowed by 
the small number statistics available. 

(ii) Evolution of the FP. At given Re and a, the in- 
termediate redshift galaxies have brighter effective surface 
brightness with respect to the values predicted by the local 
relation (solid line) and hence are more luminous and have 
lower AI/L. The average offset increases with redshift, as 
shown in panels (d). 

To make these statements more quantitative, let us de- 
scribe the evolution of the intercept 7 (at fixed slopes a and 
/3) with a linear relation: 

A7 = Tz. (22) 

With a least fit (see Figure |^ applied to the data in 
the V band we obtain r = 0.44 ± 0.03 (classical cosmology) 
and T = 0.56 ±0.04 (A cosmology); in the B band we obtain 
r = 0.45 ± 0.04 (classical cosmology) and r = 0.57 ± 0.04 
(A cosmology). This description allows us to estimate the 
scatter of the FP in any redshift bin. In fact, there are three 
contributions to the scatter: the measurement error, the in- 
trinsic scatter, and the evolutionary drift. The last item is 
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Table 1. Summary of observed quantities. For each galaxy (gal) we list redshift, central velocity dispersion (it, if measured; in km s~^), 
effective radius (re; in arcsec), and effective surface brightness {SBe; in mag arcsec"^) in the rest frame B and V bands. Details on 
the measurements are given in PIL Note that the errors on the effective radius and the effective surface brightness are correlated: the 
combination that enters the FP, re — PSBe, is better determined than the individual quantities (for a discussion see PlI; see also Figures ^ 
and 
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Figure 2. As in Figure in the A cosmology. 

simply the combined efltect of the width of the redshift bins 
and the evolution of the intercept. We can correct for the 
scatter induced by the evolutionary drift by 'evolving' the 
galaxies to the average redshift of the bin with Equation |22[ 
The corrected scatter thus obtained is remarkably small and 
constant with redshift (0.08-0.09 in 7). In Figure ^ we plot 



Figure 3. As in Figure |]j (classical cosmology), data in the rest- 
frame V band. The FP of Coma, drawn as a solid line, is taken 
from Luccy ct al. (1991). 



the evolution of the intercept 7 with redshift together with 
the best-fitting line. The thin error bars represent the total 
scatter, while the thick error bars represent the scatter cor- 
rected for evolution. It is harder to separate the contribution 
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Figure 5. Evolution of the field FP in the B (upper panels) and V 
bands (lower panels), in the classical cosmology. The symbols are 
as in Fig ^. The Fundamental Plane is seen here in the projection 
most natural for the description of errors, with photometric and 
kinematic measurements on different axes. The errors on the pho- 
tometric parameters are correlated; the combination that enters 
the FP is particularly robust (see PII). 



from measurement errors, because they are correlated and 
include a systematic component. However, by subtracting in 
quadrature the average uncertainty on 7 (for a single galaxy 
this is 0.05-0.06) one can estimate the intrinsic scatter to 
be approximately 0.05-0.07. Given the limited number of 
galaxies available per redshift bin, the error on the scatter 
is quite large (~ 30 per cent for a Gaussian distribution, 
neglecting systematics) . Nevertheless, a substantial increase 



Figure 7. Evolution of the intercept 7 of the FP and its scatter 
with redshift. The solid line is a least linear fit to the data. 
The thin error bars are the measured scatter, the thick error bars 
are t he s catter corrected for its evolutionary component (see Sec- 
tion O]). The classical cosmology is used in panels (a) and (c), 
the A cosmology in panels (b) and (d). 



in the scatter can be ruled out. Larger samples are needed 
in order to measure the scatter as a function of redshift ac- 
curately. 



4.2 The SBe-Re relation at intermediate redshift 

We will no w s tudy the evolution of the SBe-Re introduced 
in Section 2.3 on the larger sample available without inter- 
nal kinematic information. The sample includes 25 galaxies 
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Figure 8. SBe-Re relation for field early-type galaxies at inter- 
mediate redshift in the V band. The solid line is the best fit to the 
Coma data by Lucey et al. (1991). The classical cosmology has 
been used to derive effective radii. The average offset is plotted 
in panel (d) as a function of redshift. Error bars are the standard 
deviation of the mean of the sample. 



- * 


2<0.25 _ 








,x , , , ; 


0.5 1 


1.5 i 


i°s R=. - 


log hso 


1 1 1 1 1 1 1 1 
• 


1 1 1 1 1 1 




0.34<2 _ 










0.5 1 


1.5 ; 




log I'so 




0.1 0.2 0.3 0.4 0.5 



Figure 9. As in Figure H, in the A cosmology. 



out to 2; = 0.49. The results in the V band are plotted in 
Figures ^ and |^. The results in the B band are very similar 
(see figures in Section ^) . 

Following the analysis of the previous section, we de- 
scribe the evolution of the intercept of the SBe-Re relation 
(a2 of Equation with a linear relation, 

Aa2 = t'z. (23) 

By means of a least fit, we find r' — —1.30 (classical 



cosmology) and r' = —1.59 (A cosmology) in the V band. 
In the B band we find —1.37 and —1.67 respectively. 

In the context of passive evolution of stellar popula- 
tions, with the assumption that the stellar mass is propor- 
tional to the effective mass A'/, oc M, Equations |l| and || 
yield: 



< Alog(Af,/L) >= 



TZ 

2^' 



(24) 



Thus, if we assume that Re and ai (Equation ^l|) are con- 
stant with redshift, and that the evolution of the SBe-Re 
relation is given by passive evolution of M^/L, i. e. 



Aaa = ASBe = 2.5 < Alog(Af./L) >, 
r' should be related to r by: 
r = -r'/3. 



(25) 



(26) 



The values found for r and r' are in excellent agreement 
with this relation. Naturally, the agreement does not prove 
that either result is correct, since the comparison depends 
on rather strong assumptions. Nevertheless, this provides an 
obvious consistency check for both the observations (that 
rely on independent local relations) and the pure passive 
evolution interpretation. 



5 SELECTION EFFECTS 

For a physical interpretation of the results it is crucial to 
have a full control of the relevant selection effects. In partic- 
ular, from the point of view of the evolution of stellar pop- 
ulations it is important to quantify the effects of the colour 
selection and the magnitude limit. As described in PII, we 
selected our galaxies from the MDS catalog^, according to 
I8i4 < 19.3 and 0.95 < V606-l8i4 < 1.7. In Figure |l^we 
plot magnitudes and colours of the galaxies in the observed 
sample, together with the selection criteria represented by 
solid lines. 

Two main potential biases can be introduced by our 
selection. On the one hand, at high redshift (2 ^ 0.4), the 
typical luminosity approaches that of the selection criterion 
in I8i4. On the other hand, at low redshift (z = 0.1 — 0.2), 
the colours of the galaxies in our sample approach the lower 
limit in the selection criterion. In this section we will focus 
mostly on the former potential bias, while the latter will 
be better addressed in the next section, with the help of 
evolutionary population synthesis models. 

Some examples clarify these selection effects. Let us first 
consider a toy population of galaxies with the following prop- 
erties. For given a and Re, the scatter in SBe results from 
measurement errors and the intrinsic scatter of the FP; the 
colour is not correlated with the offset of the intercept. For 
this toy population, the magnitude limit will tend to select 
the brightest objects, thus simulating a stronger evolution of 
the intercept of the FP. The magnitude of the bias depends 
on the measurement error and the intrinsic scatter of the FP. 



* By comparing our photometry with that of the MDS group we 
find an average difference of less than 0.01 both in colour and Igi4 
magnitude and an r.m.s. scatter of 0.11 and 0.12 respectively in 
colour and I814 magnitude. The scatter is taken into account in 
the treatment of the selection effects. 
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Figure 10. Left panel: observed Igi4 magnitude, as measured 
by the MDS group, as a function of redsliift. Tlie limit for sam- 
ple selection (I814 < 19.3) is plotted as horizontal line. Right 
panel: observed V606-l8i4 colour, as measured by the MDS 
group, as a function of redshift. The limits for sample selection 
(0.95 <V606— 1814< 1-7) are plotted as horizontal lines. 



This effect can be quantified as follows. Let us assume tfiat 
at 2 = 0.4 (the highest redshift, where the effect is stronger) 
field early-type galaxies define an FP with slopes a and /3 
equal to those of Coma. The built-in intercept is offset by 
A7th with respect to the FP of Coma, and the scatter is 
kept as a parameter. A population of galaxies is created by 
a routine implementing a Montecarlo algorithm, the magni- 
tude limit f8i4 < 19.3 applied, and the observed offset A^oba 
recovered. This is repeated for a range of values of A'yth- 

In Figure ^ panel (a), we plot the ratio A'jobs/ Ajth 
as a function of A'yots and for various values of the intrinsic 
scatter. As expected, the bias is important for low values 
of A7 and for high values of the scatter. For our observed 
values {A^obs ~ 0.2 — 0.25, scatter ~ 0.08) the effect is small 
5 per cent ). In addition, we can use the same simulation 
to check whether the selection effects introduce any bias in 
the determination of the scatter. In panel (b) we plot the 
ratio of the recovered scatter (rmsots) to the input scatter 
(rmsth) as a function of A^fobs- For Ajots ~ 0.2 — 0.25 the 
systematic effect is much smaller (< 10 per cent) than the 
random error. 

In a more realistic second example, let us consider a 
population of early-type galaxies with fixed metallicity and 
a broad distribution of ages. At 2: ~ 0.4 our magnitude limit 
would tend to select the brightest (i.e. the youngest and 
bluest) objects. At z = 0.1 — 0.2, the colour selection crite- 
rion will pick the reddest (i.e. the oldest and faintest) ob- 
jects. In this case the colour and magnitude selection effects 
are connected, and some modeling of the evolution of the 
galactic spectra is required to study the bias. The next sec- 
tion reports a quantitative analysis of these problems based 
on a Bayesian-Montecarlo approach. 



6 CONSTRAINTS ON THE EVOLUTION OF 
STELLAR POPULATIONS 

In this section we analyse the results presented above in 
terms of passive evolution of the stellar populations. First, 
in 3.1, we compare the data to the prediction of single- 
burst passive evolution models based on Bruzual & Chariot 
(1993; GISSEL96 version, hereafter BC96) synthetic spec- 
tra. In this first qualitative analysis selection biases are ne- 



Figure 11. Estimate of the magnitude selection bias. Panel (a): 
due to the magnitude limit (see Fig. px| ) the value of A7 recov- 
ered from observation {A^y^bs) can be higher than the input value 
(A'yth)', here we plot the ratio of the observed to input value as 
a function of observed value. Panel (b): The magnitude limit can 
also bias the observed scatter (rmSoja); here we plot the ratio of 
the observed scatter to the input scatter (rmsth)- The curves are 
derived with Montecarlo simulations at z = 0.4, i.e. the highest 
redshift bin, where the effects are stronger. The curves are la- 
belled with the intrinsic scatter in 7 of the FP. The A cosmology 
is assumed. Similar results are obtained for the classical cosmol- 
ogy- 



glected. In the following Section 6.2 we introduce a Bayesian- 



Montecarlo approach (see below and the description of the 
algorithm in Treu 2001), useful to derive quantitative re- 
sults, taking into account the selection process. In its sim- 
plest version, the approach is used to calculate the a poste- 
riori probability density (i.e. the probability density given 
the present set of observations) of the t parameter defined 
in Equation Then the Bayesian-Montecarlo approach is 
combined with the BC96 spectral synthesis models to derive 
constraints on the star formation history of the galaxies in 
the sample. 

As is well known, the problem of determining the star 
formation history of stellar populations from the Fundamen- 
tal Plane (and the SBe-Re relation) is sensitive to the value 
of the cosmological parameters and the modeling of the stel- 
lar populations. In particular, it depends on the choice of 
the Initial Mass Function (IMF; see Schade et al. 1997). For 
this reason, the evolution of the FP can be used to con- 
strain the cosmological parameters and/or the IMF (Bender 
et al. 1998; van Dokkum et al. 1998b). In the present study, 
in order to limit this sort of degeneracies, we focus on the 
'best-guess' for the value of the cosmological parameters and 
the form of the IMF (A cosmology and Salpeter IMF) and 
then we try to extract the maximal information on the star 
formation history. In order to give an idea of what is robust 
and what is sensitive to the assumptions about the cosmo- 
logical parameters and the IMF, we will also present the 
results for the classical cosmology and for t he Scalo IMF. 
For completeness we briefiy discuss (Section 6.3) the limits 
that can be set on the cosmological parameters. 

In Sec tion 3.4 we discuss potential systematic biases and 



in Section 6.5 we extend th e d iscussion of the population 
biases introduced in Section 2.2, 
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Figure 12. Single-burst passive evolution models. The evolution 
of the intercept of the FP is compared to the prediction of single- 
burst passive evolution models computed from BC96 synthetic 
spectra in panels (a) and (b). The comparison is repeated in pan- 
els (c) and (d) for the SBe-Re relation, by including galaxies 
without measured velocity dispersion. The solid lines represent 
models with Salpeter IMF (Salpeter 1955), the dashed lines rep- 
resent models with Scalo IMF (Scalo 1986). All models have solar 
metallicity. Three redshifts of formation are assumed, Zf = 1,2,5, 
from top to bottom in panels (a) and (b), from bottom to top in 
panels (c) and (d). The classical cosmology is assumed. 



6.1 Single-burst models 

In Figures and |l| the evolution of 7 and 02 as pre- 
dicted by single-burst passive evolution models is shown 
and compared to the data. For each cosmological model we 
show the predictions obtained for single bursts occurred at 
Zf = 1,2,5. The results are shown for models with solar 
metallicity and kl96 stellar atmospheres (see the documen- 
tation to BC96 for details). The results are highly insensitive 
to changes in metallicities and stellar atmospheres. 

Qualitatively, we note that the data points are better 
reproduced by models with Zf Jj 2 in the classical cosmology 
and with 2;/ ~ 1 in the A cosmology. 



6.2 A Bayesian-Montecarlo approach to the 
evolution of stellar populations 

We can better quantify the results obtained so far by means 
of a Bayesian-Montecarlo approach. In the following we will 
give only a brief description of the method, which is outlined 
in detail elsewhere Treu (2001). 



6.2.1 The evolution of the intercept of the FP 

Let us assume that EquatioriO holds and that the inter- 
cept 7 evolves as in Equation p2[i.e. linearly in redshift with 
slope T (A7 = Tz). Under these assumptions we are able to 
calculate the probability — p*(7'|t) of observing a galaxy 
at redshift with 7' (defined as in Equation [l^), where the 



Figure 13. As in Figure |l2| in the A cosmology. 

superscript i indicates an individual galaxy within the set of 
observed galaxies. The probability density is obtained from 
Montecarlo simulations (at least 1,000,000 random events 
per galaxy are generated, see Treu 2001), taking into ac- 
count the intrinsic scatter of the FP (assumed to be constant 
with redshift as measured above), the observational errors, 
and the selection effects. Since galaxies are independent, the 
probability of observing the set {7'} is: 



P({7'}k) = (7». 



(27) 



We now apply Bayes Theorem (Bayes 1763; see also, e.g.. 
Zellner 1971): 



Pi-r\{l'}) = 



P{W}\r)p{r) 
J p{{l'}\'r)p{T)dT' 



(28) 



to derive the probability of a value of r given the set of 
measured {7'} (the a posteriori probability, or simply the 
posterior). The probability density of r prior to observa- 
tions (the a priori probability, or simply the prior) can be 
used to take into account previous observations or phys- 
ical limits. The choice of the a priori probability density 
in cases where investigators know little about the value of 
the parameter, or wish to proceed as if they knew little, 
is extensively discussed in the literature (see, e. g., Zellner 
1971 and references therein; see also the recent review by 
Berger 1999). Traditionally, in the so-called Bayes-Laplace 
approach, a uniform prior was used. Alternative priors have 
been proposed, including some based on the requirement of 
invariance with respect to changes of parametrization (Jef- 
freys 1961) and others based on information arguments (see 
next paragraph; see also Zellner 1971). The advantage of 
the Bayesian approach is that it provides a probability dis- 
tribution for the parameter (as compared for example to the 
confidence limits provided by a 'frequentist' approach), un- 
der well-controlled assumptions (such as the choice of the 
prior) . 

For simplicity, in the following we assume p{t) constant 
within the interval (—0.1, 1.5) and zero outside it, i.e. 
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Figure 14. Evolution of tlie intercept of tlie FP in the V band. 
The a posteriori probability density of the slope r of the relation 
A7 = Tz is shown here as recovered from a Bayesian-Montecarlo 
approach that takes into account selection biases (see Section 3.2). 



The posteriors obtained with a uniform prior (solid line) and with 
a prior of the form p(t) oc 1/r (dashed line) are shown. The results 
for the classical cosmology and the A cosmology are shown in 
panels (a) and (b). 



p{t) = 0.625J/(1.5 - r)H{r + 0.1), 



(29) 



with H the Heaviside function, defined as H{x) — 1 for 
3; > and otherwise. The interval has been chosen so as to 
include all the region where p({7'}|t) 7^ 0. It is interesting to 
note that for a finite interval spanned by the parameter, the 
uniform prior is the proper (i.e. normalizable to 1) prior with 
minimal information (defined as J p{x)\ogp{x)dx; Zellner 
1971). It is also noticed that for this choice of the prior the 
maximum of the posterior is the value that would be derived 
with a maximum likelihood algorithm. 

The method described above, extended to the galaxies 
without measured velocity dispersion in terms of the SBe- 
Re relation as outlined by Treu (2001), gives the results 
shown in Figure |l| (V band). One sigma limits on r can be 
defined in the usual way, by requiring that the integral of 
the probability density over the interval be 0.68. In this way, 
we find 0.33 < t < 0.44 (classical cosmology) and 0.44 < 
r < 0.56 (A cosmology). These intervals are consistent with 
the values found using a least fit (see Equation |22| and 
the following discussion). However, the values found with 
the least fit are on the higher end of the interval, not 
on the peak of the probability density (see Figure We 
conclude that a correct treatment of the selection effects is 
important in order to derive reliable estimates. 

In order to check that the results are robust with respect 
to the choice of the prior, we computed the posterior by 
assuming an alternative prior. For example, let us consider 
p(t) oc 1/t (limited to r > 0), which has the advantage 
that the results are unchanged under a power law change 
of parametrization X = r" (see Zellner 1971 and references 
therein). In fact, by converting 



p{X\{Y})^pi{Y}\X)p{X)<x 
into a posterior in r we get 

p'{r\W})<xp{Xir)\W}) 



X 



dX 



dr 



Pi{Y}\X{r)) 



(30) 



(31) 



as is obtained by assuming p(r) cc 1/r. The posterior ob- 
tained by taking this prior is shown in Figure M as a dashed 



line; it is very similar to the one obtained by assuming a uni- 
form prior (solid line). 

6.2.2 Single-burst models 

The Bayesian-Montecarlo approach is used with the BC96 
modeling of integrated stellar populations to obtain quanti- 
tative information on the star formation history of our sam- 
ple of field early- type galaxies. Let us first consider a single 
population of galaxies with fixed IMF and stellar popula- 
tions formed in a single-burst at z = 2/ and evolving pas- 
sively thereafter. Let us assume that in the local Universe 
such population obeys the FP relation of Coma. As usual, 
we will assume that the properties of the stellar populations 
are the only factor inducing evolution. Therefore, only the 
intercept of the FP evolves with redshift, while the scatter 
and slopes a, /3 remain constant. In this model, for any given 
redshift, velocity dispersion, efi'ective radius (2*, cr*, rg), and 
observational errors of an individual galaxy, we can calcu- 
late the probability density of the colours, rnagnitude, and 
intercept of the FP (7' defined in Equation |l2| ) . By taking 
into account the selection effects with a Montecarlo routine 
(a large number of points are generated and only the ones 
satisfying the selection criteria are kept) we obtain the prob- 
ability density p'(7'|2/) of observing 7* given Zf. Since the 
galaxies are independent samplings of the statistical ensem- 
ble we obtain: 



p{H}\zf) = ny(7>/). 



(32) 



As in Equation the probability density of Zf given the 
set of observations can be computed by assuming an a priori 
probability density p{zf) from the Bayes theorem: 



P{zf\{j}) = 



J Pi{'y'}\zf)pizf)dzf' 



(33) 



In the following, we will assume a constant a priori proba- 
bility between two values of 2: 



p{zf) = plH{zf — Zm)H{zM — Zf). 



(34) 



A very natural choice is one with Zm very close to the highest 
2' of the sample and a high value of zm, in order to include 
the interval where p(7*[2/) is non-zero. We use Zm ~ 0.5 and 
Zm = 4.0 for the A cosmology and Zm ~ 0.5 and zm = 10.5 
for the classical cosmology. 

The a posteriori probability densities are shown in Fig- 
ures |l|and0 The probability density peaks at 2 ~ 1 — 1.2 
(8-9 Gyr look-back time) and it extends to 2 ~ 0.8 (7 Gyr) 
and 2 ~ 1.8 (11 Gyr) in the A cosmology. In the classi- 
cal cosmology the probability is much flatter and it extends 
from 2 ~ 1.5 — 2 (10 Gyr) to the upper limit of the sam- 
pled interval (the age of the Universe in this cosmology is 
13.1 Gyr). The distributions are computed for a variety of 
realizations of the model, with different metallicities, IMFs, 
stellar atmospheres, showing that the result is robust and 
does not depend on the details of the model. 

In order to check the dependence of the results on the 
a priori probability density we consider the following ap- 
proach. An alternative and equivalently straightforward a 
priori probability is the one obtained by considering the 
look-back time of formation (tf) as the parameter of the 
model, instead of the redshift of formation (zf) used in 
Equation pa: 
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Figure 15. Single-burst stellar population models, I. The proba- 
bility density of the redshift of formation (zf) of the stellar pop- 
ulations of field early-type galaxies, given the present set of ob- 
servations, is shown for various realizations of the models. BC96 
spectral synthesis models are used with Salpeter or Scalo IMF, so- 
lar or 0.4 solar metallicity, stellar atmospheres from kl96 or from 
gsHR (see Bruzual & Chariot 1996 for details). Independently of 
the model, the probability density turns out to peak at Zf ~ 1. 
The A cosmology is assumed. Details of the method are given in 
Section 3.2 



p{tf\W}) 



piW}\tf)pitf) 
J p{{r}\tf)p{tf)dtf' 



(35) 



Similarly to the case with Zf, we assume a prior uniform 
within the interval defined by the look-back time of the most 
distant galaxy in our sample and the age of the Universe. 
The posterior p(t/|{7*}) can be transformed into the poste- 
rior in in the usual way: 



p'(^/|{7})=p(t/(2/)|{7}) 



dtf 



dzf 



where 
dtf 



dzf 
with 



Ho{l + Zf)y/n{l + Zf)3 + Qr(1 + ZfY + Qa 



(36) 



(37) 



Smce p({7%^) = p{{l'}\tf{zf)), 
tions B3, B5l and ha it follows that: 



P'iZflW}) «P(2/1{7'}) 



dtf 



dzf 



(38) 

from Equa- 
(39) 



i.e. the probability density obtained by assuming a prior 
uniform in t/ is proportional to the one obtained by assum- 
ing a prior uniform in z/ multiplied by the modulus of the 
derivative of tf with respect to z/ (the constants are set by 
requiring that the probability density integrated over Zf he 
equal to 1). 



Figure 16. Single-burst stellar population models, II. As in Fig- 
ure for the classical cosmology. 



The resulting a posteriori probability densities are plot- 
ted in Figures |l7| and respectively for the A and the 
classical cosmology. It is noticed that the a posteriori prob- 
ability density changes negligibly in the A cosmology, while 
in the classical cosmology the probability density peaks at 
smaller Zf and declines at high Zf. However, this does not 
alter the physical interpretation of the result, since in any 
case the stellar populations turn out to be formed long ago 
(^ 10 Gyr). Since the effects of changing the a priori proba- 
bility density do not alter the physical interpretation of the 
results, in the following we will present only the results ob- 
tained by assuming a prior uniform in z/, recalling that the 
results obtained by assuming a prior uniform in tf can be 
obtained with Equation 

6.2.3 Secondary bursts of star formation 

The single-burst stellar population model is useful as a 
benchmark to quantify and compare different results. How- 
ever, it is clearly a simplified picture of the star formation 
history of early- type galaxies. A small amount (< 10 per 
cent) of the total mass of moderately young stars (1 Gyr or 
so) can significantly alter the integrated colours and mass- 
to-light ratio of an old stellar population for a few Gyr. In 
this scenario, after a few Gyr, the integrated colours and 
mass-to-light ratio become totally indistinguishable from 
the ones of an old single-burst stellar population (see also 
Jimenez et al. 1999). 

A wealth of observations suggest that minor episodes of 
star formation can occur from intermediate redshifts to the 
present (e.g. Schade et al. 1999; Trager et al. 2000a; Bernardi 
et al. 1998) . Therefore, it is worth investigating this possibil- 
ity for our sample of field early-type galaxies. The task can 
be accomplished by modifying the model described above 
to allow for a secondary burst of star formation occurred at 
Zf2 < 2/1, being Zfi associated with the epoch of the first 
burst. In this case Equation is changed to 
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Figure 17. Single-burst stellar population models, III. Effects 
of the a priori probability density. The a posteriori probability 
density of the look-back time of formation (tf) and of the red- 
shift of formation (zf) are plotted respectively in the upper and 
lower panel (A cosmology). Dashed lines correspond to a prior 
uniform in tj-, solid lines to a prior uniform in zj^. BC96 spectral 
synthesis models are used with Salpeter IMF, solar metallicity, 
and kl96 stellar atmospheres. The effects of changing the a priori 
probability are negligible. A di scussion of the a priori probability 
density is given in Section 6.2. 




1 + z, 



Figure 18. Single-burst stellar population mod< 
ure hTj, for the classical cosmology (see Section 



la^IV. As in Fig 
3.2). 



Figure 19. Model with two populations of stars. The older com- 
ponent is formed at Zfi, the younger, with a tenth of the mass, 
is formed at Zf2- Contour levels of the a posteriori probability 
are shown as solid (68 per cent) and dashed (95 per cent) lines. 
The A cosmology is assumed. BC96 models with solar metallicity, 
Salpeter IMF and kl96 atmospheres are used. 



P(^/l.^/2|{7'}) = 



Pi{l''}\Zfl,Zf2)p{Zfl,Zf2) 
J P{{T}\zfl,Zf2)p{zfl, Zf2)dZfidZf2 



.(40) 



The a priori probability is assumed to be constant within a 
test interval 

p{zfl,Zf2) OC H{zfl — Zm)H{zM - Zfi) X 

H{Zf2 - Zm)H{zM - Zf2)H{Zfl ~ Zf2), (41) 

similarly to the one-dimensional case. The contour levels of 
the probability density corresponding to 68 per cent and 
95 per cent probability are shown in Figure 1^ for the A 
cosmology and a mass ratio of the two populations of 10 
(the older being 10 times the mass of the younger). Not 
surprisingly, it is sufficient to have a small mass of stars 
formed at z < 0.6 — 0.8 to make it possible for the rest of 
the stellar mass to be formed at the very beginning of the 
Universe {zfi ~ 3 — 4, i.e. 12-13 Gyr ago). 

It is important to check that such scenario is also con- 
sistent with the small scatter of the FP observed at low and 
intermediate redshift. Since a number of physical phenom- 
ena contributes to the scatter of the FP, the scatter due to 
differences in ages must be smaller than the observed one. 
In order to simulate this effect we have considered a very 
simple model, where the bulk of the stellar population is 
formed at z — zt and a secondary burst (of a tenth of the 
mass) occurs at a later time, with uniform probability be- 
tween z = Zsi and z — Zs2- The observed trend in A7 vs. z 
is reproduced without any particular fine tuning of the pa- 
rameters {zb ~ 2 — 3, Zsi ~ 0.6, Zs2 ~ 0.8, for example, work 
well) and the resulting scatter in 7 is very small (^ 0.01). It 
is interesting to note that, if we extend to lower redshifts the 
tail of the secondary burst of star formation, the evolution 
of 7 is not altered significantly, while the scatter induced 
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by the spread in average age of the stellar populations is 
comparable to the observed one 0.08). 



6.3 Constraining the cosmological parameters 

The Bayesian-Montecarlo approach can be extended also to 
constrain the cosmological parameters. To this aim, it is 
sufficient to modify Equation into 



p{{^'}\zf,n)p{zf,n) 
J p{{r}\zf,n)p{zf,n)dzfdn' 



(42) 



where, for simplicity, only Q has been introduced, and to 
assume a discrete a priori probability using the Dirac 5 dis- 
tribution 

p{zf, Q) = piH{zf - Zml)H{zMl - Zf)S{Q - 0.3) 

+P2H{zf - Z,^2)H{ZM2 - Zf)5{Q - 1). (43) 

The normalization factors (pi, P2) are chosen so that the in- 
tegral of each term over Zf and SI equals 1/2. The intervals 
in Zf are chosen as in Section 3.2. The probability density 
pdy^Mzf, O.), is related to the probability densities in Equa- 
tion by 



p{m\zf,n) oc pn=o.3i{Y}\zf)S{n - 0.3) 
+pn=i{W}Uf)m^i)- 



(44) 



The probabilities obtained for the two cosmologies are com- 
pared in Figure The maximum of probability occurs 
clearly for the A cosmology but the area subtended by the 
classical cosmology is actually larger (0.3 vs 0.7). Hence, at 
this redshift, given the present uncertainties on the evolution 
of the stellar populations, the A cosmology is the 'best-guess' 
cosmology, but the classical cosmology cannot be ruled out 
by this method (see the 68 per cent limits marked with thick 
lines in Figure EOt see also the discussion by Bender et al. 
1998 and van Dokkum et al. 1998b). 



N 



T3 




5 6 7 



Figure 20. Probability density of the redshift of star-formation 
given the present set of observations. The a posteriori probabil- 
ity densities for the classical and the A cosmology are compared 
by means of an extension of the Bayesian-Montecarlo approach 



(see Section 5.3). The maximum of the probability is for the A 
cosmology. However, the 68 per cent limits (thick lines) include 
also the classical cosmology. 



of Coma and the FP of other clusters or the field. We also ne- 
glect possible small differences between the FP in the centre 
and in the outer parts of Coma (Lucey et al. 1991). 



6.4 Systematic errors 

The analysis performed above takes the random errors of 
the measurement into account. However, we still have to 
estimate the uncertainty related to global systematic errors, 
such as: 

(i) Zero point errors. We estimate the uncertainty in 
the absolute zero point of filters F606W and F814W to be <, 
0.03 mag (see the HST Instrument Handbook and references 
therein for the listing and description of uncertainties) 

(ii) Errors in the K-color correction. These errors 
are ~ 0.01 — 0.02 mag (see PH). As a check, we compared 
the K-color corrections obtained with our method with the 
ones derived by Kelson et al. (2000a) and van Dokkum et 
al. (1998b). We found, consistent with the error estimate, 
differences ~ 0.01 — 0.02 mag. 

(iii) Uncertainty in the local FP. The slopes and the 
zero point of the FP of the Coma cluster are known with er- 
rors. The uncertainty on 7 for our reference samples (Bender 
1998, Lucey 1991) is ~ 0.01, including photometric calibra- 
tion, K-correction, and foreground extinction. We assume 
the distance to Coma to be cz = 7200 km s~^ and we ne- 
glect, as discussed in Section 0, differences between the FP 



The total systematic error, obtained by adding in 
quadrature the terms listed above, does not exceed 0.015 
in 7. We tested the robustness of our results with respect to 
this uncertainty by recomputing the probability density for 
T after changing the Coma intercept by ±0.015. The 68 per 
cent limits are changed by ±10 per cent. 

Since the coefficients of the FP and SBe-Re relation 
are measured independently, the very good agreement of the 
results found in B and V, with the FP and SBe-Re relation, 
is an important check on the accuracy of the photometric 
calibration and the determination of the local relations. 

In principle, it is possible to take into account system- 
atic errors in the Bayesian-Montecarlo approach. In fact, it 
is sufficient to treat the photometric zero point, the K-color 
corrections, and the intercept of the local FP as additional 
parameters of the model and then integrate over the values 
of the parameters, as in the following example. Let us de- 
note the set of additional parameters related to the system- 
atic errors as and let us assume that their probabil- 
ity distribution are independent p{{Y^}) = nip(y). Thus, 
Equation ^ becomes: 

^ J piH}\T{Y'})pir)U,pHY')dY' 
P[mi i) Jp(^{Y}\T{Y'})p{r)n,p^{Y')dY'dT' ^ ' 



© 0000 RAS, MNRAS 000, 000-000 



16 T. Treu et al. 



However, the computational cost of implementing this pro- 
cedure in our Bayesian-Montecarlo code|T] is not justified in 
our opinion, given the limited knowledge of the systematics, 
and the other sources of uncertainties that enter the inter- 
pretation of the results (see also the next subsection) . 



6.5 Morphological evolution bias 

The morphological selection criterion can introduce a bias 
in our res ults (the morphological evolution bias described 
in Section 2^). For example, one may think (Franx & van 
Dokkum 1996) that early-type galaxies become morphologi- 
cally recognizable only after their stellar populations reach a 
minimum age (Aq). In this case, the evolution of the morpho- 
logically selected sample would be biased towards an older 
average age of the stellar population (simply because objects 
younger than Aq are discarded). 

A simple model where E/SO form at a random time in 
a given redshift interval can help quantify the bias. For ex- 
ample, in the A Universe, by forming galaxies in the redshift 
range z = 0.5 — 1.5, and by considering them as E/SO only 
after 3 Gyr, we can reproduce the observed evolution of 7 
with constant scatter ~ 0.06 — 0.07. 

However, in this scenario no E/SO would be recognizable 
at 2 = 1, in contrast with the observations (Im et al. 1996; 
Schade et al. 1999; Treu & StiaveUi 1999). By requiring that 
at leEist one third of present-day E/SO be already assembled 
and recognizable at z = 1, we may need a smaller Ao or 
an earlier period of formation. Reducing ^0 does not work, 
since the age component in the scatter of the FP rapidly 
exceeds the values allowed by the observations. An earlier 
period of formation is therefore needed. Starting earlier the 
formation (z ~ 3 is required) alone does not work either, 
since the evolution of 7 becomes too slow. In order not to 
exceed the observed scatter it is necessary to stop forming 
E/SO earlier in the Universe. For example, if the formation 
occurs between 2 = 1 and 2 = 2 and Aq — 2 Gyr, the 
evolution of 7 is reproduced, the age contribution to the 
scatter is always smaller than the scatter observed, and 48 
per cent of E/SO are recognizable at 2 = 1. 

In conclusion, based on the results of the FP alone, 
we cannot exclude that field E/SO have been forming until 
recently (2 = 0.5) if they become morphologically recogniz- 
able only at a moderately old age (~ 3 Gyr). In turn, we can 
exclude this possibility, if we take into account the observa- 
tional constraint that a significant part of the population of 
E/SO is recognizable at 2 = 1. With this constraint we are 
forced to push back the epoch of formation of the stellar 
populations of field E/SO to 2 ^ 1. 



t In practice the CPU time (a few hours on a Pentium III running 
at 500 Mhz for the computations shown in Figures ^ to [l7[ ) is 
multiplied by the number of samplings of the distribution of the 
parameters related to systematic errors. Even by taking only ten 
samplings per systematic parameter (V), the CPU time increases 
1000 times. 



Table 2. Compilation of sources for FP parameters of galaxies 
in clusters at intermediate redshift. For each cluster we list the 
redshift, the rest frame band in which photometry is given, the 
number of objects per cluster, and the reference. Only objects 
classified as E/SO have been considered. See Section [t] for discus- 
sion. 
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Figure 21. Comparison of field and cluster ellipticals in the FP- 
space at intermediate redshift. Filled pentagons are field data 
points from this paper, onen squares are cluster data taken from 
the literature (see Table bh. In panels (a) and (c) the position in 
the FP-space of galaxies at similar redshift is shown. In panels 
(b) and (d) the offset of the intercept with respect to the value of 
Coma is plotted as a function of redshift. The A cosmology is as- 
sumed. The comparison yields very similar results in the classical 
cosmology (see Treu 2001). 



7 COMPARISON BETWEEN FIELD AND 
CLUSTER RESULTS 

In this section we compare the intermediate redshift field FP 
to a compilation of cluster data taken from the literature (see 
Table 1). 

In Figure ^ we plot the cluster and the field FP data 
points in the A cosmology (the results of the comparison are 
very similar for the "classical" cosmology and are not shown 
here, see Treu 2001). In panels (a) and (c) the cluster and 
field data points at 0.33 < z < 0.42, respectively in the V 
and B band, are plotted in the FP-space. Panels (b) and (d) 
illustrate the offset of the intercept with respect to Coma as 
a function of redshift. The same slopes a and /3 have been 
used for consistency, thus eliminating this source of system- 
atic error in the relative comparison. The data at similar 
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1+z, 

Figure 22. Comparison of tlie redshift of formation probability 
density obtained for field (solid line) and cluster E/SO (dashed 
line). The cluster data in the B band have been considered (see 
Table |2|). BC96 models with Salpeter IMF, solar metallicity, and 
kl96 atmospheres have been used. The A cosmology is assumed. 
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Figure 23. As in Figure Only cluster galaxies in the same 
ranges of a and Re as those of the field sample have been used 
(open squares). The results found with the complete cluster sam- 
ple are plotted as crosses for reference. 



redshift {z ~ 0.3 — 0.4) provide the following results. The 
clusters points at 2 = 0.375 and z = 0.39 (panels b and d) 
are marginally (1-2 standard deviations) below the field data 
points, while the point at z = 0.33 is well consistent with 
the field data points. By applying the Bayesian-Montecarlo 
code to recover the a posteriori probability density for the 
cluster data points, we confirm that: the stellar populations 
of the galaxies in the two clusters ai z — 0.375 and z — 0.39 
are marginally older than the ones in the field. On the other 
hand, the probability distribution obtained for the galaxies 
in the cluster at z = 0.33 is indistinguishable from the one 
obtained in the field. We conclude that, given the differences 
in the selection process (see below), and the size of observa- 
tional errors, no difference is detected between cluster and 
field galaxies. 

The cluster data points at higher redshift {z — 0.58 — 
0.83) show consistently a smaller offset from the local rela- 
tion than the one that would be extrapolated from the field 
points. For example, if we apply our Bayesian-Montecarlo 
approach to the data in B band out to z = 0.83, we obtain 
an older average age for the cluster stellar populations (see 
Figure . A larger sample of field data, extended to higher 
redshift, is needed in order to perform a more stringent com- 
parison. To this aim, we are currently collecting spectra for 
25 galaxies in the redshift range z = 0.4 — 0.7 using FORS2 
on the Very Large Telescope (VLT). We expect that the ob- 
servational errors on 7* will be comparable with the ones 
obtained for the present study. Thus, we will be sensitive to 
possible differences of ~ 10 per cent between the offset of 
the intercept 7 for field and cluster E/SO. 



7.1 Sample selection effects 

The galaxies in the cluster sample span wider ranges in a 
and Re than those of our sample of galaxies. This effect is 
related to the lower limit in velocity dispersion measurable 
by our instrumental setup (see PIT) and by the absence of 
giant objects such as brightest cluster galaxies in our sample. 
Note that we were not biased against very large (and hence 
bright) objects. We simply did not find any of them. This 
selection effect could bias the results if, for example, larger 
galaxies were older than smaller ones. 

The size of this bias can be estimated by comparing 
the field data to a subsample of cluster E/SO with a and 
7?e within the limits of our sample. This subsample shows a 
marginally larger offset from the local relation than the com- 
plete sample (see Figure ^) . However, this result should be 
taken with great care, since the samples we are considering 
are very small and the change in the intercepts can be simply 
due to noise. Larger samples with well controlled selection 
effects are needed to address the issue whether the evolution 
of the stellar populations depend on effective mass, possibly 
with the help of a measurement of a change in the slopes of 
the FP with redshift. 

In addition, the possible bias induced by the cluster 
population evolution (see Section ^ has to be considered. 
Phenomena such as the Butcher-Oemler effect (Butcher & 
Oemler 1984), the evolution of the morphology-density re- 
lation (Dressier et al. 1997), or the larger scatter of the 
colour-magnitude relation of lenticular galaxies, noticed by 
van Dokkum et al. (1998a) in the outer parts of a cluster 
at z = 0.33, seem to suggest that the galaxy population of 
clusters evolves. Therefore, it is difficult to extend the infor- 
mation gathered by studying the FP in the core of clusters 
at high redshift to the whole population of present-day clus- 
ters. In particular, if galaxies are accreted from the field into 
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the cluster, the cluster population at intermediate redshift 
is not the progenitor of the present-day cluster population. 
Studies with a large field of view, encompassing the entire 
range of densities from the core of rich clusters to the sparse 
random fields, arc needed before wo can clarify how cluster 
populations evolve and how the environment affects galaxy 
evolution. 



8 CONCLUSIONS 

In this paper we have studied the Fundamental Plane and 
the SBe-Re relation for a sample of field eaxly-type galaxies 
at intermediate redshift. In particular, we have used the FP 
and SBc-Rc as diagnostics of the stellar populations. The 
main results are: 

(i) Evolution of the FP. The FP of field early-type 
galaxies exists and is tight out to « w 0.4. The scatter is 
consistent with being unchanged out to z ^ 0.4. The inter- 
cept is offset with respect to that of Coma, in the sense that 
at a given effective radius and velocity dispersion galaxies 
arc brighter at z ~ 0.4 than in the local Universe. The ranges 
of a, SBe, and re covered by our sample are not sufficient to 
measure the slopes of the FP at intermediate redshift accu- 
rately. Similar results are obtained by studying the SBe-Re 
relation out to z ~ 0.5. 

(ii) Single-burst stellar populations. In a single-burst 
scenario, the observed properties of this sample of galaxies 
are consistent with those of a single-burst stellar population 
formed at z ~ 0.8 — 1.6 (A cosmology) or at z J> 2 (classical 
cosmology). These redshifts of formation correspond to a 
present age of 7-11 Gyr or 10-13 Gyr, respectively, in the 
two cosmologies. 

(iii) Multiple-burst stellar populations. If a small 
fraction of the stellar mass is formed in a secondary burst at 
a later time, the primary burst may have occurred at high z 
even in the A cosmology. For example, the data are consis- 
tent with the scenario where the primary burst occurred at 
z ^ 3 and a secondary burst, with a tenth of stellar mass, 
occurred at 2: ~ 0.6 — 0.8. 

(iv) Field vs. Cluster. No significant difference is found 
between the Fundamental Plane of field E/SO at intermedi- 
ate redshift and that of cluster E/SO at similar redshift taken 
from the literature. The ages of the stellar populations of 
field early-type galaxies inferred from the single-burst model 
are marginally smaller than the ages derived for samples 
of cluster early- type galaxies extended to higher redshift. 
Unfortunately, various selection processes, the small sam- 
ples available, and the different ranges in global properties 
spanned by the field and cluster samples do not allow us to 
draw firm conclusions. Larger sets of data at higher redshift, 
as the one we are currently collecting with the VLT, are 
needed to investigate possible differences with higher sensi- 
tivity. A larger set of data, sampling a larger range of effec- 
tive masses and a variety of environments from rich clusters 
to the field, is also needed to address the issue of the clus- 
ter population evolution bias and the possible dependence 
of the evolution on effective mass. 
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